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This viewpoint article is intended as a brief introduc- 
tion to the emerging subject of field-theoretic simulations 
(FTS) of charged polymers. While the direct numerical 
simulation of field theory models has begun to impact 
several traditional areas of polymer science, including 
blends and block copolymers, polyelectrolytes have hith- 
erto not been the subject of field-theoretic simulations. 
Here we report on a preliminary FTS study of polyelec- 
trolyte complexation that demonstrates the potential of 
this novel numerical approach. 

Polyelectrolytes are ubiquitous in nature and in appli- 
cations ranging from personal care products to paints, 
coatings, and processed foods. Indeed, practically all 
biopolymers are polyelectrolytes. In the application con- 
text, the introduction of dissociable groups is one of 
the most powerful ways to confer water solubility on a 
polymeric material. Scientifically, the polymer bound 
charges, which are compensated by a sea of oppositely 
charged counterions, produce a coupling between chain 
conformations and electrostatics that leads to an incredi- 
ble richness of polyelectrolyte phenomena. However, this 
richness comes with a price: charged polymers are among 
the most difficult polymer systems to study theoretically 
or to simulate on the computer [l|, 0, H, Q . The expla- 
nation lies in the long range nature of Coulomb interac- 
tions - charged segments feel each other at much larger 
distances than segments in neutral polymer systems. 

The situation becomes particularly challenging for 
the dense polyelectrolyte complexes that are the sub- 
ject of the present work. Analytical theories based on 
assumptions of low concentration or weak interactions 
break down, and equilibration times in numerical simu- 
lations become prohibitively long. To conduct such sim- 
ulations, polyelectrolytes are often modeled as coarse- 
grained chains of charged beads and the counterions are 
taken to be point particles, usually embedded in an im- 
plicit solvent with uniform dielectric properties. The 
long-ranged character of Coulomb interactions is prob- 
lematic for such "particle-based" modeling approaches, 
however, because Ewald sums and other expensive com- 
putational techniques are required to evaluate contribu- 
tions to electrostatic energies and forces that extend be- 
yond the computational cell d, H, The confluence 
of this long-range effect with the intrinsically slow ki- 
netics of particle-based simulations at high density and 
molecular weight produces significant challenges for these 
methods. 
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A different strategy to tackle such problems was pro- 
posed long ago by S. F. Edwards [8]. The idea is to re- 
place the coordinates and momenta of particles (polymer 
segments) with collective variables, or fields. One can 
introduce, for example, fields that describe the density 
of polymer segments and the density of charge. More- 
over, by augmenting these fields with certain conjugate 
fields, such as a chemical potential (conjugate to seg- 
ment density) and an electrostatic potential (conjugate 
to charge density), it is possible to exactly transform any 
classical (equilibrium) particle-based model into a sta- 
tistical field theory This field-based description is 
particularly useful for dense polymer systems, such as 
concentrated solutions and melts, where there is strong 
overlap among polymers. In such cases, the mean-field 
approximation, also known as self-consistent field theory 
(SCFT), can be applied with confidence. Within this ap- 
proximation, the saddle point of the Hamiltonian dom- 
inates the partition function, and the field fiuctuations 
around the saddle point are ignored The SCFT 

method has yielded some remarkable analytical and nu- 
merical results for a wide variety of dense equilibrium 
polymer systems, perhaps most notably in the field of 
block copolymers [ll|, [l^ . 

Besides its restriction to systems at equilibrium, a ma- 
jor limitation of SCFT is that the assumption of negli- 
gible field fluctuations breaks down rapidly as polymers 
are diluted in a solvent. This is especially problematic for 
the study of polyelectrolytes, since they are highly sol- 
vated in most applications. Moreover, polyelectrolytes 
are characterized by very strong charge correlations in 
addition to the density correlations of neutral polymer 
solutions. Evidently the well-tuned machinery of SCFT 
is not the appropriate tool for investigating solutions of 
charged macromolecules. 

One way to address the limitations of SCFT, while still 
preserving the field-based collective variable approach, is 
to return to the exact statistical field theory and to ac- 
count for the field fluctuation effects. Analytically, this 
procedure leads to a so-called "loop" expansion, in which 
systematic corrections to the free energy or other ther- 
modynamic quantities (one-loop, two-loop, etc.) are de- 
veloped in terms of integrals over field fluctuations about 
the saddle point of the theory. Such loop expansions can 
be further augmented by renormalization techniques in 
cases of strong fiuctuations 1, 13]. While powerful in the 
context of homogeneous phases of charged and uncharged 
polymers, these analytical methods can be very difficult 
to implement in more general inhomogeneous situations 
where interfaces or mesophases are present. 

It has recently been demonstrated that statistical field 



theory models of polymers can also serve as the basis for 
computer simulations - so-called field-theoretic simula- 
tions (FTS) [1, While similar conceptually to nu- 
merical SCFT, field-theoretic simulations aim to numer- 
ically sample the statistically important field configura- 
tions of the full theory, rather than just the saddle point 
configuration. This statistical sampling is problematic 
because the Hamiltonians of the relevant field theories 
are complex, rather than strictly real; a manifestation 
of the famous "sign problem". We have found that a 
powerful way to circumvent this difficulty is to adopt the 
complex Langevin stochastic procedure (isl . [l6| . which 
adaptively samples field configurations along nearly con- 
stant phase trajectories. Although FTS is computation- 
ally more expensive than SCFT, recent numerical ad- 
vances have made high resolution FTS feasible for a wide 
variety of polymer systems [l7| . 

Polyelectrolytes have been the subject of extensive the- 



oretical and computational research for decades 
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22, [S^l ■ Statistical field theory models have played a sig- 
nificant role in these theoretical investigations, and both 
mean-field and non-mean-field approaches have been em- 
ployed to gain insights into the structure and ther- 
modynamics of a wide variety of polyelectrolyte sys- 
tems HH, m, m [13, m m, H . Prior to this work, 
however, there has been no general numerical tool for 
simulating a field theory model of polyelectrolytes with- 
out the use of simplifying approximations. 

The present article describes the first application 
of FTS to polyelectrolytes, and specifically to a phe- 
nomenon that is known to occur in aqueous mixtures 
of two oppositely charged polymers. Under appropriate 
conditions of charge density, solvent quality, and molec- 
ular weight, a phase transition can occur in such a sys- 
tem, with two phases being formed: one rich in both 
polymers and the other consisting of nearly pure sol- 
vent [H, [H [13, [H, [H, [13, [M [si, [13. This process 
is usually referred to as 'polyelectrolyte complexation in 
the physics community or as complex coacervation in the 
physical chemistry, colloid science, and biological com- 
munities. The resulting polymer-rich coacervate phase 
has two important properties: it is dense yet liquid, and 
it is charge neutral. 

Coacervates and related polyelectrolyte complexes 
have a variety of important applications. For example, 
complexes can be employed as carrier systems for charged 
macromolecules including protein drugs, enzymes, and 
DNA. In such systems, one of the polyelectrolytes serves 
as a chaperon with properties that can be tailored to 
assist targeted delivery, while the other represents the 
macromolecular payload Another industrial use of 
complex coacervation is in cements, glues, or adhesives, 
where the coacervate phase can serve as a precursor to 
the formation of a solidified structural material. Ex- 
amples of such coacervate precursors are also encoun- 
tered in nature; e.g. sand-castle worms produce a strong 
protein-based cement that sets in sea water and is used 
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FIG. 1: Two model polyelectrolyte mixtures. Both are sym- 
metric, with half of all chains being positively charged (shown 
in black) and the other half negatively charged (in grey), (a) 
The total charge ±aN is evenly distributed over the entire 
length of each chain, (b) The same total charge is concen- 
trated in blocks of N/A: segments near one end of each chain. 
Uncharged portions are shown in white. 



to construct meter-scale reefs out of sand and shells [41| . 
Other existing and potential applications of polyelec- 
trolyte complexes include water purification, DNA sen- 
sors [1^, and encapsulation of food and pharmaceuticals. 

Polyelectrolyte complexation is driven by several com- 
peting factors [li,[li,l33,[3i,[3i,[33,[3i,[3i|. In addition 
to the direct electrostatic attraction between oppositely 
charged polyions, there are other factors playing an im- 
portant role: electrostatic screening by small ions (salt 
or counterions) , excluded volume interactions of polymer 
backbones mediated by the solvent, small ion transla- 
tional entropy, and polymer conformational entropy. In 
many cases these factors are competitive rather than re- 
inforcing, so it can be difficult to anticipate the conditions 
under which coacervate phases will form. For example, 
excluded volume interactions oppose complexation, while 
the translational entropy of counterions tends to favor it. 

To investigate the physics of complex coacervation 
we have chosen a simple yet fundamental model sys- 
tem [Fig. [ija)]. Specifically, we consider a symmetric 
polycation-polyanion mixture in an implicit solvent with- 
out salt. Such a system could be realized by mixing a 
polyacid with a polybase in water. The symmetric as- 
sumption implies that the molecular weights and charge 
densities of the two types of polymers are the same, and 
they are combined in equal amounts. The polymer back- 
bones are modeled as fiexible (Gaussian) chains of length 
(number of statistical segments) N , differing only by the 
sign of the charge per statistical segment (charge den- 
sity), ±cr. We employ a canonical ensemble in which n 
polyanions and n polycations are mixed to form a solu- 
tion of volume V . The chains are assumed to interact by 
means of electrostatic and excluded volume interactions, 
characterized by the Bjerrum length Ib = e^/eksT and 
the excluded volume parameter uq, respectively. The sol- 
vent dielectric constant is denoted by e, and e is the fun- 
damental unit of charge. Note that the addition of salt, 
consideration of asymmetric polyions, or explicit inclu- 
sion of the solvent constitute only minor modifications 
of the formalism, although these variations will not be 
pursued here. 

The thermodynamic properties of the above model can 
be deduced from a configurational partition function Z, 



3 



which is a sum over all possible states (represented by 
the coordinates of all polymer segments) of the Boltz- 
mann factor exp(— [//fcsT). The energy of each state U 
can be expressed as a sum of three contributions: Gaus- 
sian chain connectivity, excluded volume interactions, 
and electrostatic (Coulomb) interactions. Rather than 
working with this "particle-based" statistical mechanical 
description, we reexpress the partition function in terms 
of collective field variables 0. Indeed, Z can be written 
as a sum (a functional integral) over all field configura- 
tions: 



Z^Zo 



J I DwD(f) exp{-H[w, 



(1) 



where Zq is the partition function of an ideal gas of non- 
interacting polymers, w{r) is a real chemical potential 
field conjugate to the segment density, and (/)(r) is a real 
electrostatic potential field conjugate to the total charge 
density. The energy of each state is now replaced by an 
effective Hamiltonian, which for our model is 
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In this effective Hamiltonian, direct interactions between 
polymers decouple, and the remaining interactions are 
field-field and field-polymer. Polymers interact with the 
fields only via the single- chain partition function Qlip], 
where "0(1") corresponds to the two purely imaginary 
fields iw ± ia(f) in EqO This partition function is calcu- 
lated as a volume integral of a chain propagator, which 
in turn can be obtained from the solution of a complex 
diffusion equation in the imaginary field ip f9] . Kj^ is the 
coefficient in front of the Green's function of the Lapla- 
cian in d dimensions (47r in 3D and 27r in 2D). 

It is important to note that this field theory is an ex- 
act reformulation of the original particle-based model. 
While the fundamental model is not new, previous 
field-theoretic approaches to polyelectrolyte complexa- 
tion [23, m, [131 have simplified the field theory by use 
of weak inhomogeneity expansions. This simplification is 
not exact, nor is it desirable or necessary. We further 
note that the long-ranged Coulomb interaction in the 
particle description is replaced by a short-ranged "square 
gradient" interaction |V0p in the field theory. From a 
computational standpoint, this is very attractive. The 
field-theoretic representation is also convenient for an- 
alytical studies of field fluctuation effects, such as loop 
expansions. 

A convenient place to begin our discussion of the field 
theory model is with the mean-field approximation, or 
equivalently, SCFT. The mean-field solution neglects all 
field fluctuations and replaces the sum over fleld con- 
flgurations in Eq [1] by a single most probable conflgu- 
ration - a saddle point located off the real axis in the 
complex plane. For the present model (subject to peri- 
odic boundary conditions) the saddle point configuration, 



obtained from the simultaneous equations SH/5w = 
and 5H/5(j) = 0, is homogeneous (a constant pure imagi- 
nary number) in both fields. Furthermore, due to global 
charge neutrality, the saddle point solution has a distinc- 
tive feature: the mean electrostatic potential is constant 
so that every positive charge is compensated by an equal 
negative charge. It follows that the Coulomb interactions 
are irrelevant in the mean-field limit and the polyelec- 
trolyte mixture should behave exactly as an analogous 
mixture of neutral polymers. This leaves open no pos- 
sibility for complexation, and indeed we shall see that 
it is necessary to include fiuctuations and account for 
charge correlations to obtain a coacervate phase. Thus, 
the mean-field approximation breaks down completely in 
the case of a simple mixture of cationic and anionic poly- 
electrolytes. 

Analytically, the next level of sophistication is to ac- 
count for quadratic field fluctuations about the saddle 
point solution in the evaluation of Eq[T] This is the lead- 
ing term in a systematic loop expansion (l3| known as 
the Gaussian or one-loop approximation. Calculations 
of this type have been reported previously for models of 
polyelectrolyte complexation [13, [H, [131 ■ For our sym- 
metric polycation-polyanion mixture, fluctuations of the 
w and (j) fields decouple at the one-loop level. Moreover, 
it turns out that the system is parameterized by only 
three dimensionless combinations of variables: reduced 
polymer concentration C = 2nRg/V, reduced excluded 
volume parameter B — uqN^ / Rg, and reduced Bjerrum 

length E = Kahcr^ / Rg~^ ■ Here Rg = {NlP /2dY/^ is 
the size of an ideal non-interacting polymer (which coin- 
cides with its radius of gyration in 3D), h is the statistical 
segment length, and d is the number of dimensions. Thus, 
the one-loop approximation provides relevant combina- 
tions of parameters responsible for the thermodynamic 
state of the system [43]. It should be noted that the 
electrostatic interactions are manifested only in the pa- 
rameter E. 

The one-loop approximation yields analytical expres- 
sions for correlation lengths in the polyelectrolyte mix- 
ture. In particular, it predicts that polymer segment 
density fluctuations are correlated on the scale of the 
well-known Edwards length ^ . Of greater interest is the 
prediction that charge density fluctuations are correlated 
on the scale of the polyelectrolyte length 
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(3) 



where p = 2nN/V is the segment density. This cor- 
relation length, which has been previously identified in 
3D [H, [13], has several important characteristics. First, 
it is specific to charged polymers, as can be observed from 
its dependence on the Bjerrum length and the statisti- 
cal segment length. Second, S^pe is independent of the 
chain length N , as is expected for dense polymer sys- 
tems. Finally, it is proportional to the —1/4 power of 
the segment density and hence is qualitatively different 



4 



from the Debye-Hiickel length ^dh = {^/ KdlBC^ pY^^ 
for smaU ions of density p carrying charge ±ct. Thus, the 
attachment of charges to polymer chains creates a cou- 
pling between chain conformational statistics and charge 
density that dramatically changes the electrostatic cor- 
relation properties of the solution compared with a con- 
ventional small ion electrolyte. 

Beyond correlation lengths, the one- loop approxima- 
tion provides an analytical expression for the Helmholtz 
free energy that can be used to derive expressions for 
standard thermodynamic quantities. For example, the 
fluctuational electrostatic contribution to the osmotic 

pressure scales as - Ad{EC)'^/^kBTR-'^ kBT/£_f,j^, 

where Ad > is a known numerical coefficient dependent 
on dimensionality d. This expression is again qualita- 
tively different from the analogous electrostatic term in 
Debye-Hiickel theory for simple electrolytes. Because the 
net contribution from charge correlations is negative, the 
one-loop theory can predict polyelectrolyte complexation 
and the coexistence of dilute and coacervate phases. The 
spinodals and binodals for this two-phase region follow 
directly from the one-loop osmotic pressure expression. 
While some of these predictions were obtained previously 
by slightly different methods [H, HI, [131 , their validity 
has been difhcult to assess because the next term in the 
loop expansion (two-loop order) is very tedious to evalu- 
ate. 

With the advent of the FTS method, we now have a 
numerical technique that can be used to simulate the full 
field theory without the assumption of weak charge and 
density correlations that underpins the loop expansion. 
Computer simulations of the full field theory can pro- 
vide test beds for the analytical loop expansions, just as 
particle simulation methods complement analytical virial 
expansions in the particle-based theory. 

The advantage (or disadvantage) of FTS over conven- 
tional particle-based simulation techniques can be as- 
sessed by comparing computational costs. The cost of 
particle-based simulations depends on the total number 
of atomistic or coarse-grained particles in the system. 
State-of-the-art particle-based methods require a num- 
ber of operations per MD step or MC cycle on the order 
of nNlnnN for n polyelectrolytes, each with TV beads or 
segments @ . On the other hand, the cost of an FTS field 
update depends on the spatial discretization of the sys- 
tem, rather than the number of particles. When the com- 
putational cell is divided into a lattice of size M, each up- 
date requires of order NM In M operations @ . In a melt, 
M < n if the lattice spacing Ax in the FTS approach 

1 /3 

can be taken larger than Vp , where ~ iV is the vol- 
ume of a polymer. This is met, for example, by choos- 
ing Ax to be a fraction of Rg - N^/'\ Thus, FTS has 
an obvious computational advantage for dense systems 
of long polymers in which a coarse computational grid 
suffices to capture mesoscopic structure at and beyond 
the Rg scale. On the other hand, particle-based methods 
are advantaged under more dilute conditions, when fewer 
molecules need to be described on microscopic scales, or 



when full chemical details are required. These considera- 
tions are rough guides, and further studies are needed to 
fully elucidate the conditions under which FTS is com- 
petitive with existing simulation techniques. 

We employ complex Langevin (CL) sampling in our 
FTS to avoid the sign problem associated with the com- 
plex Hamiltonian H and the non-positive definite char- 
acter of the statistical weight exp(— i?) in the field the- 
ory. This method was originally devised as a strategy for 
sampling general types of quantum field theories with 
complex actions ^15i . ,44] and has been more recently ap- 
plied in polymer physics [ll,|4l,|4^. The idea behind this 
method is to extend the real fields into the complex plane 
and to compute ensemble averages of observable quanti- 
ties by sampling fields along a stationary stochastic tra- 
jectory in the complex function space. While this exten- 
sion to the complex plane doubles the number of field de- 
grees of freedom, the relevant statistical weight becomes 
a real non- negative distribution of fields P[W, $], where 
W = wji + iwi and $ = 0^ -I- i^/. 

A stochastic CL dynamics in the complex func- 
tion space is employed to generate a Markov se- 
quence of complex fields with stationary distribution 
P\W, $]. The complex Langevin equations are dW/dt — 
—A ((5iJ/(5M^) -I- 77 and a similar expression for the $ field. 
Here W and 811 jSW (or <i> and 5H/5^) are complex, but 
the thermal noise 77(r, t) is a real Gaussian white noise, 
satisfying the usual fluctuation-dissipation theorem with 
dissipative coefficient A. Since the thermal noise is placed 
asymmetrically only on the real part of the CL dynamics, 
the imaginary parts of the equations ensure the sampled 
fields have nearly constant phases Hj along the Langevin 
trajectory. This removes rapid oscillations and improves 
convergence. Recent improvements in stochastic integra- 
tion algorithms for the CL equations have enabled high- 
resolution, three-dimensional (3D) field-theoretic simula- 
tions 

Here we report on the application of CL-FTS to poly- 
electrolyte complexation phenomena, specifically to the 
same symmetric binary polycation-polyanion model that 
was used for the analytical one- loop calculations. A par- 
ticular focus of our CL-FTS study is the location and size 
of the two-phase region where coacervation takes place. 
An example of such a phase diagram (in 3D) is provided 
in Fig. O Spinodals and binodals for the field theory of 
Eq[2]are surfaces in the three-parameter space of the re- 
duced variables C, B, and E. The figure represents a 
cross-section of this three-dimensional space by a plane 
E = 14400; hence the diagram involves only the C and B 
variables. The diagram features both the CL-FTS results 
(symbols and dotted-line fit) and the one-loop analytical 
approximation to the binodal (solid line) and the spin- 
odal (dashed line). The one-phase region (disordered ho- 
mogeneous phase) is above and to the right of the lines, 
and the two-phase region is below and to the left. The 
tie lines in the two-phase region are horizontal (constant 
B) and connect nearly pure solvent (C ~ 0) with a coac- 
ervate phase at the binodal concentration. The two sym- 
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FIG. 2: Phase diagram for the symmetric polyelectrolyte mix- 
ture in 3D plotted in reduced polymer concentration C and 
reduced excluded volume B at fixed reduced Bjerrum length 
E. Solid and dashed lines are the analytical one-loop binodal 
and spinodal, respectively. Symbols are the result of com- 
plex Langevin simulations; the dotted line is a power- law fit. 
Numerical simulations were conducted in a cubic cell of size 
47?g X ARcj X ARg with periodic boundary conditions. 



bols for the CL simulation data at each concentration C 
represent the hysteresis upon varying the solvent quality 
B, with the upper symbol corresponding to superheat- 
ing (increasing B) and the lower one corresponding to 
supercooling (decreasing B). 

Remarkably, the analytical and numerical results 
nearly coincide in the high concentration region of the 
figure, despite the limitations of each method. The nu- 
merical results are subject to finite cell size and chain dis- 
cretization limitations, while the analytical predictions 
neglect two-loop and higher order terms in fiuctuations. 
Nonetheless, our numerical supercooling result practi- 
cally follows the analytical spinodal, and the analytical 
binodal yields nearly the same exponent (—1.31) as is ob- 
tained from a power-law fit to the numerical superheating 
points (—1.40). The cause for the discrepancy between 
theory and simulation for the location of the binodal is 
unclear at present, although the overall semi-quantitative 
agreement indicates that both approaches have utility for 
this class of problems. 

Our future work in polyelectrolyte complexation will 
extend beyond the symmetric model presented here to 
include unequal chain lengths, counterions, salt, and ex- 
plicit solvent. The latter will allow for the treatment 
of polyelectrolytes with hydrophobic backbones. It is 
important to emphasize that these are straightforward 
extensions that do not complicate the particle-to-field 
transformation of the model, nor do they make analytical 
loop expansions or CL simulations any more difficult. 

Another extension is to investigate the effect of charge 
distribution and polymer architecture on polyelectrolyte 
complexation. This brings us into the realm of block 
copolyelectrolytes and branched and dendritic polyelec- 
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(a) 




(b) 

FIG. 3: Snapshots of equilibrated field-theoretic simulations. 
Normalized total polymer density is presented for the cases of 
(a) uniformly charged chains and (b) unevenly charged chains, 
with the same amount of total charge concentrated on 25% of 
the chain length. Macrophases are formed in the uniformly 
charged case, and a mesoscopic structure develops in the case 
of block polyelectrolytes. Both systems are in 2D with peri- 
odic boundary conditions. Cell size is I6-R3 x l&Rg, and the 
same parameters C — 6.0, B — 0.3, and E — 64000 were used 
in both cases. 



trolyte systems - all of considerable experimental inter- 
est. As a preliminary example, we have explored a sim- 
ple variation of our binary polycation-polyanion model 
where all the charge of each species is concentrated into 
a small "charged block" at the chain end; c.f. Fig. [TJb). 
The remaining longer portion of each chain constitutes a 
"neutral block" . In our simple model, the charged blocks 
should drive complexation because of the correlation- 
induced electrostatic attraction, while the neutral blocks 
repel each other due to the excluded volume interactions 
in the (assumed) good solvent. These competing ten- 
dencies can drive aggregation to form micellar structures 
with charged blocks in the micelle cores and the neutral 
blocks in the coronas. At relatively low concentrations, 
the micelles are expected to form a disordered micellar 
phase, while at higher concentrations they can pack into 
periodic lattices to form mesophases of various symme- 
try. Such structured coacervate phases would seem to be 
of considerable interest for a variety of applications. 
In accord with expectations, complex Langevin simu- 
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lations of the symmetric binary block copolymer model 
show that the redistribution of charges along the poly- 
electrolyte chains leads to the formation of aggregates 
qualitatively similar to micelles. Figure [3] illustrates 
the dramatic difference between the phases formed by 
uniformly charged chains and unevenly charged chains 
with the same total charge concentrated into a block of 
length N/A. The uniformly charged chains [Fig. El^a)] 
form macrophases, i.e. large homogeneous phases (a di- 
lute phase and a coacervate phase) constrained only by 
the size of the computational cell. On the other hand, 
the oppositely charged block copolymers [Fig. ^h)] form 
a mesophase with structure on the scale of Rg. The core 
of each micelle predominantly consists of charged blocks, 
while the neutral blocks form the repulsive corona. We 
are aware of at least one body of experimental work that 
confirms these (as yet) qualitative predictions [47l. [isj. 

In summary, field-theoretic simulation methods based 
on the introduction of auxiliary fields have considerable 
promise for the investigation of polyelectrolyte complexa- 
tion phenomena. These methods possess several distinc- 
tive and attractive features. First, FTS involves exact 
Hamiltonians and accounts for arbitrarily large fluctua- 
tions and strong inhomogeneities. Thus the technique is 
not limited in the same ways as the available analytical 
tools. Second, field-theoretic methods are especially con- 
venient for treating the long-range Coulomb interaction, 
which is replaced by a short-ranged square-gradient op- 
erator in the auxiliary field representation. Third, FTS 



methods are computationally advantaged over particle- 
based simulations when the systems are dense, polymer 
chains are long, and the length scale of interest is large 
(mesoscopic) . Finally, inclusion of additional species (ei- 
ther long chains or small ions) does not lead to a substan- 
tial elaboration of the field-theoretic models, nor does it 
complicate the simulations. 

With continued advances in algorithms, we antici- 
pate that FTS methods will enable the study of broad 
classes of polyelectrolyte systems and phenomena that 
are beyond the reach of current analytical methods and 
particle-based computer simulations. The importance of 
analytical theory and particle-based simulations will not 
be diminished, however, by the emergence of FTS tech- 
niques. We see these approaches as complementary, with 
each contributing valuable insights into the rich structure 
and thermodynamics of charged polymeric fluids. 
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